Molecular Mechanism of Methanol Inhibition in CALB-Catalyzed Alcoholysis: Analyzing Molecular Dynamics Simulations by a Markov State Model.

Lipases are widely used enzymes that catalyze hydrolysis and alcoholysis of fatty acid esters. At high concentrations of small alcohols such as methanol or ethanol, many lipases are inhibited by the substrate. The molecular basis of the inhibition of Candida antarctica lipase B (CALB) by methanol was investigated by unbiased molecular dynamics (MD) simulations, and the substrate binding kinetics was analyzed by Markov state models (MSMs). The modeled fluxes of productive methanol binding at concentrations between 50 mM and 5.5 M were in good agreement with the experimental activity profile of CALB, with a peak at 300 mM. The kinetic and structural analysis uncovered the molecular basis of CALB inhibition. Beyond 300 mM, the kinetic bottleneck results from crowding of methanol in the substrate access channel, which is caused by the gradual formation of methanol patches close to Leu140 (helix α5), Leu278, and Ile285 (helix α10) at a distance of 4-5 Å from the active site. Our findings demonstrate the usefulness of unbiased MD simulations to study enzyme-substrate interactions at realistic substrate concentrations and the feasibility of scale-bridging by an MSM analysis to derive kinetic information.

Discovery of phosphotyrosine-binding oligopeptides with supramolecular target selectivity.

We demonstrate phage-display screening on self-assembled ligands that enables the identification of oligopeptides that selectively bind dynamic supramolecular targets over their unassembled counterparts. The concept is demonstrated through panning of a phage-display oligopeptide library against supramolecular tyrosine-phosphate ligands using 9-fluorenylmethoxycarbonyl-phenylalanine-tyrosine-phosphate (Fmoc-FpY) micellar aggregates as targets. The 14 selected peptides showed no sequence consensus but were enriched in cationic and proline residues. The lead peptide, KVYFSIPWRVPM-NH2 (P7) was found to bind to the Fmoc-FpY ligand exclusively in its self-assembled state with K D = 74 ± 3 µM. Circular dichroism, NMR and molecular dynamics simulations revealed that the peptide interacts with Fmoc-FpY through the KVYF terminus and this binding event disrupts the assembled structure. In absence of the target micellar aggregate, P7 was further found to dynamically alternate between multiple conformations, with a preferred hairpin-like conformation that was shown to contribute to supramolecular ligand binding. Three identified phages presented appreciable binding, and two showed to catalyze the hydrolysis of a model para-nitro phenol phosphate substrate, with P7 demonstrating conformation-dependent activity with a modest k cat/K M = 4 ± 0.3 × 10-4 M-1 s-1.

Integration of Molecular Dynamics Based Predictions into the Optimization of De Novo Protein Designs: Limitations and Benefits.

Recent advances in de novo protein design have gained considerable insight from the intrinsic dynamics of proteins, based on the integration of molecular dynamics simulations protocols on the state-of-the-art de novo protein design protocols used nowadays. With this protocol we illustrate how to set up and run a molecular dynamics simulation followed by a functional protein dynamics analysis. New users will be introduced to some useful open-source computational tools, including the GROMACS molecular dynamics simulation software package and ProDy for protein structural dynamics analysis.

Comparison of the Internal Dynamics of Metalloproteases Provides New Insights on Their Function and Evolution.

Metalloproteases have evolved in a vast number of biological systems, being one of the most diverse types of proteases and presenting a wide range of folds and catalytic metal ions. Given the increasing understanding of protein internal dynamics and its role in enzyme function, we are interested in assessing how the structural heterogeneity of metalloproteases translates into their dynamics. Therefore, the dynamical profile of the clan MA type protein thermolysin, derived from an Elastic Network Model of protein structure, was evaluated against those obtained from a set of experimental structures and molecular dynamics simulation trajectories. A close correspondence was obtained between modes derived from the coarse-grained model and the subspace of functionally-relevant motions observed experimentally, the later being shown to be encoded in the internal dynamics of the protein. This prompted the use of dynamics-based comparison methods that employ such coarse-grained models in a representative set of clan members, allowing for its quantitative description in terms of structural and dynamical variability. Although members show structural similarity, they nonetheless present distinct dynamical profiles, with no apparent correlation between structural and dynamical relatedness. However, previously unnoticed dynamical similarity was found between the relevant members Carboxypeptidase Pfu, Leishmanolysin, and Botulinum Neurotoxin Type A, despite sharing no structural similarity. Inspection of the respective alignments shows that dynamical similarity has a functional basis, namely the need for maintaining proper intermolecular interactions with the respective substrates. These results suggest that distinct selective pressure mechanisms act on metalloproteases at structural and dynamical levels through the course of their evolution. This work shows how new insights on metalloprotease function and evolution can be assessed with comparison schemes that incorporate information on protein dynamics. The integration of these newly developed tools, if applied to other protein families, can lead to more accurate and descriptive protein classification systems.

A bis(pyridine N-oxide) analogue of DOTA: relaxometric properties of the Gd(III) complex and efficient sensitization of visible and NIR-emitting lanthanide(III) cations including Pr(III) and Ho(III).

We report the synthesis of a cyclen-based ligand (4,10-bis[(1-oxidopyridin-2-yl)methyl]-1,4,7,10-tetraazacyclododecane-1,7-diacetic acid=L1) containing two acetate and two 2-methylpyridine N-oxide arms anchored on the nitrogen atoms of the cyclen platform, which has been designed for stable complexation of lanthanide(III) ions in aqueous solution. Relaxometric studies suggest that the thermodynamic stability and kinetic inertness of the Gd(III) complex may be sufficient for biological applications. A detailed structural study of the complexes by (1) H NMR spectroscopy and DFT calculations indicates that they adopt an anti-Δ(λλλλ) conformation in aqueous solution, that is, an anti-square antiprismatic (anti-SAP) isomeric form, as demonstrated by analysis of the (1) H NMR paramagnetic shifts induced by Yb(III) . The water-exchange rate of the Gd(III) complex is ${k{{298\hfill \atop {\rm ex}\hfill}}}$=6.7×10(6)  s(-1) , about a quarter of that for the mono-oxidopyridine analogue, but still about 50 % higher than the ${k{{298\hfill \atop {\rm ex}\hfill}}}$ of GdDOTA (DOTA=1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid). The 2-methylpyridine N-oxide chromophores can be used to sensitize a wide range of Ln(III) ions emitting in both the visible (Eu(III) and Tb(III) ) and NIR (Pr(III) , Nd(III) , Ho(III) , Yb(III) ) spectral regions. The emission quantum yield determined for the Yb(III) complex (${Q{{{\rm L}\hfill \atop {\rm Yb}\hfill}}}$=7.3(1)×10(-3) ) is among the highest ever reported for complexes of this metal ion in aqueous solution. The sensitization ability of the ligand, together with the spectroscopic and relaxometric properties of its complexes, constitute a useful step forward on the way to efficient dual probes for optical imaging (OI) and MRI.